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1 Introduction 



A wide range of cosmological observations indicate that the universe has entered a 
phase of accelerating expansion, which becomes one of the important puzzles of the 
contemporary physics. Those observations include the type la supernovae (Snla) 
standard candles PQ, [2J, the angular location of the first peak in the CMB power 
spectrum [3] and baryon acoustic oscillations of the matter density power spectrum 
[I] . These evidences represent a great stimulus for theoretical work and originated the 
concept of Dark Energy. The Dark Energy (DE) models attribute the observed ac- 
celerating expansion to an unknown energy component with negative pressure, which 
dominates the universe at recent cosmological times. The simplest model describing 
the cosmic expansion rate with a high degree of accuracy, is the cosmological con- 
stant [SJ, [U] in a spatially flat space, which provides a constant equation of state. 
However this physically motivated model is plagued by fine tuning problems [7j. Al- 
ternatively, dark energy may be described by different dynamical scalar field models 
with time-dependent equation of state, like quintessence [8], [9]; string theory funda- 
mental scalar known as tachyon [TU]; K-essence models involving a generalized form 
of the kinetic energy [H],[T2]; scalar field with negative kinetic term, which provides 
a solution known as phantom dark energy [13] (see [H] for a review). An alternative 
description of DE may be given by perfect fluids with adequate equation of state, 
like Chaplygin gas [15J. Another possibilities to the explanation of the DE are rep- 
resented by the scalar-tensor theories, which contain a direct coupling of the scalar 
field to the curvature, providing in principle a mechanism to evade the coincidence 
problem, and naturally allowing (in some cases) the crossing of the phantom barrier 
|lbj . The non-minimal coupling between the quintessence field and curvature, as an 
explanation of the accelerated expansion, has been considered among others, in refs. 
[TTj [18],[19j,[20j. An inflationary model including the DBI term and non-minimally 
coupled scalar has been proposed in [2TJ [22] . 

A model that can be regarded as falling within the scalar-tensor theories is the model 
with non-minimally coupled kinetic terms to the curvature. In this paper we consider 
an explicit coupling between the scalar field, the kinetic term and the curvature, as a 
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source of dark energy, and analyze the role of this coupling in an evolution scenario 
with late-time accelerated expansion [231 121] • The basic motivation for studying such 
theories is related with the fact that they appear as low energy limit of several higher 
dimensional theories, e.g. superstring theory [22], and provide a possible approach to 
quantum gravity from a perturbative point of view [26]. A coupling between curva- 
ture and kinetic terms also appears as part of the Weyl anomaly in N = 4 conformal 
supergravity [271 128] . A model with non-minimal derivative couplings was proposed 
in [22], [30], [31] in the context of inflationary cosmology, and recently, non-minimal 
derivative coupling of the Higgs field was considered in [32] , also as inflationary model. 
In [33] a derivative coupling to Ricci tensor has been considered to study cosmological 
restrictions on the coupling parameter, and the role of this coupling during inflation. 
Some asymptotical solutions for a non-minimal kinetic coupling to scalar and Ricci 
curvatures were found in [31], and quintessence and phantom cosmological scenarios 
with non-minimal derivative coupling have been studied in [33]. A scalar field with 
kinetic term coupled to a product of Einstein tensors has been considered in |36j , and 
a vector field non- minimally coupled to curvature has been studied in [37] . A kinetic 
coupling of scalar field to Ricci curvature appears as part of the action in a class of 
covariant renormalizable gravity theories [38J. 

The modified gravity theories have received increased attention in the last years in 
connection with the dark energy problem, specifically the f(R) theories [39], [4*0] . 
0U> S3> US], n US, n H3> n n |S3- The modified Gauss Bonnet model 
has also been considered as gravitational alternative for dark energy in [21], [52] . 
A modified F(R) Horava-Lifshitz gravity and, a unified description of inflation and 
late-time acceleration in the context of modified F(R) Horava-Lifshitz gravity, have 
been considered in [23], [21], [22]- 

Non-minimal coupling of arbitrary function f(R) with matter Lagrangian (in- 
cluding the kinetic scalar term) has been introduced in [2E],[27]. Such a model was 
proposed to describe the dark energy and late-time universe acceleration. When f(R) 
represents the power law function, such a model (which maybe considered as string- 
inspired theory) maybe proposed for dynamical resolution of cosmological constant 
problem as it has been suggested in [2H],[29]- In the present work we consider this 
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function f(R) as linear in R but we generalize the model permitting extra cou- 
pling with kinetic-like scalar term. In addition, we keep the free kinetic term for the 
scalar field. 

In the present paper we continue the study of the model [21], in which the scalar-field 
kinetic term is non-minimally coupled to itself and to the scalar and Ricci curvatures, 
and in addition we consider the matter contribution. The self coupling function F ((/)), 
allows to find a wide variety of potentials giving rise to accelerated expansion. The 
power-law solution is considered, and the expressions for the scalar-field and poten- 
tial are found. The solution corresponding to the cosmological constant with matter 
(ACDM) is obtained as an exact solution in the present model, giving in this way 
a dynamical interpretation to the cosmological constant. In addition to the ACDM, 
other relevant solutions have been studied. In the next section we review the main 
equations in a general background and in the flat FRW metric. In section 3 we study 
power-law solutions giving rise to accelerated expansion. In section 4, the ACDM and 
other relevant solutions for the dark energy are studied, first in the case of scalar field 
dominance, and then including the matter contribution. Solutions that may accom- 
plish the current observations have been considered. In the last section we present 
some conclusions. 



2 Field Equations 

Let us start with the following action, proposed in [24] : 

S = J d A x^~g [j^R ~ \d^4> - ±£R (F(0)<V<9>) - ^ ^ 

^R^ (F(<f>)d»<f>d»(f>) - V + S m . 

where S m is the matter action described by a fluid with barotropic equation of state 
(along the text, if it is not specified the term "matter" will enclose the ordinary 
barionic matter and dark matter, which obey the same equation of state). The 
dimensionality of the coupling constants £ and r\ depends on the type of function 



F{4>). Taking the variation of action 2.1 with respect to the metric, we obtain a 
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general expression of the form 

R 



1 



g, u R = K 2 [T^ + T, u ] 



(2.2) 



where n 2 = 8nG, T™ is the usual energy-momentum tensor for the ordinary and 
matter component, and the tensor T^ u represents the variation of the terms which 
depend on the scalar field (f> and can be written as 



T — 4- T"£ 4- T"? 



(2.3) 



where T$ v , T^ u) TJJ U correspond to the variations of the minimally coupled terms, the 
£ and the r\ couplings, respectively with respect to the metric. Due to the interacting 
terms these expressions are defined in the Jordan frame and do not correspond to 
the energy-momentum tensors as defined in the Einstein frame. Therefore, due to 
this interaction between the scalar field and the curvature, the derived expressions for 
the density and pressure for the scalar field can be regarded as effective ones. Those 
variations are given by 



It = V^V^ - -^V A 0V A - 9llv V{<t>) 



T% H R,u - ^9,uR (F(0) V A 0V A 0) + ^,V A V A (F(0) V 7 0V» 



;(V,V„ + V,V„) (F(0)V A 0V A 0) + R (F(0)V^V„0) 



(2.4) 



(2.5) 



1% =V 



F(<j>) (i^ A V A 0V,0 + iUV>V M 0) - -g^R Xj (F(0)V>V 7 0) 



" 2 (^(0)V A 0V,0) + V A V, (F(0)V>V M 0)) 

+ ^V A V A (F(0)V^V„0) + ^V A V 7 (F(0)V A 0V^0) 
Varying with respect to the scalar field gives rise to the equation of motion 



(2.6) 



dF 

d(f> 



(2.7) 



{iRd^d^ + V R^d^d»<P) = 



Assuming the spatially-flat Friedmann-Robertson- Walker (FRW) metric, 



(2.8) 
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from Eqs. (2.4 2.6) and using the metric (2.8) we can write the (00) and (11) com- 



ponents of the Eq. (2.2) (with the Hubble parameter H, and for homogeneous time- 
depending scalar field) as follows 



H 2 



'Ptot 



and 

with p tot given by 

Ptot 



2H - 3H 2 



K 2 ptoU 



(2.9) 



(2.10) 



1 ■ 



2 + V(<t>) + 9££PF(0)0 2 + 3(2£ + n)HF 



- 3(2£ + r))HF 



(2.11) 



+ Pr, 



and ptot given by 
-1 



Ptot 



- V((f>) + 3(f + ri)H 2 F(0)4> 2 + 2(f + r))HF 



+ 4(f + v)HF 

+ (2Z + V ) (f(0)0 2 + F(0)0 



2^ + V )H^ 3 



(2.12) 



5dF 



+ 



ld 2 F • 



2 c/0 r r ' 2 d0 2 " 

here we have assumed that the matter contribution corresponds to presureless perfect 



fluid (i.e. p m = 0). The equation of motion for the scalar field (2.7) takes the form 



• dV ... . . 

+ 3H<f> +-T- + 3(2£ + rj)HF(<p)<p + 3(14f + 5r/)MF(0)0 



+ -(2£ + T?)tf(2F(0)0 + 



dF ,, 



+ -(4£ + 7 / )if 2 (2F(0)<H 



(2.13) 



+ 9(4£ + r])H 3 F(<f))<p = 
where "dot" denotes derivative with respect to cosmic time. The first three terms in 



Eq. (2.13) correspond to the minimally coupled field. In what follows we study the 
cosmological consequences of these equations under some conditions that simplify the 
search for solutions, but nevertheless show the richness of the cosmological dynamics 
of the present model. 
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First note that the Eqs. (2.9 2.13) significantly simplify under the restriction on £ 
and 7] given by 

7] + 2£ = (2.14) 
This restriction is equivalent to a coupling of the kinetic term to the Einstein tensor 



G^y (see [30], [31]), as can be seen from expression (2.1). In this case the field 
equations (2.9 2.13) contain only second derivatives of the metric and the scalar field. 



The modified Friedmann equations (|2.9|) and (|2.10|) take the form 
H 2 



K 



-0 2 + ^(0) + 9^F(0)0 2 + p, 



(2.15) 



and 

-2H-3H 2 



1 , 



(j) 2 - V{<f>) — i 3.ff + 2H ) F 



2£H 2F(0)0 4> + 



dF ;. 



(2.16) 



and the equation of motion reduces to 



<j> + 3H<P + + 3£H 2 [ 2F 
dtp 



+ 



dF 



+ 6£H (3H 2 + 2H )F(0)0 = O (2.17) 



In what follows we will study cosmological solutions to Eqs. (2.15) and (2.17), giving 
rise to accelerated expansion and that satisfy the demands of the actual observational 
data. 



3 Power- law solutions 



Let us consider the model (2.1) in the case of scalar field dominance, without matter 
contribution. Here we study solutions for the coupling function F(<p) and the poten- 
tial V(<p), which give rise to power-law expansion. 

A. Power-law solutions without potential V(4>) = 



7 



First we consider the model without potential, and analyze the role of the coupling 



function F((f>) in the power-law expansion scenario. From Eq. (2.15) it follows 

1 X2 



3£k 2 18£# 2 



(3-i; 



Making V = in Eq. (2.17), multiplying by (f> and defining the new variable <fi 2 = if}, 



the Eq. (|2.17|) reduces to first order differential equation with respect to both variable 

(3.2) 



if) and H. 







considering the power-law behavior H = s/t, the Eq. (3.2) becomes 



t 3d ^ + (Qs + l)t 2 iP + (3s-2] 



6s z 







(3.3) 



which has the general solution 

m = 



6s 2 (3s - 2) 7 



(3.4) 



where 7 is the integration constant with the appropriate dimension. This equation 
gives the following solution for the scalar field 



4>{t) 



6s 



l-e 



A + v A arcsin 



^3a-l/2 



A 



6s 2 (3s - 2) 
(6s — 



(3.5) 



and the time dependence of the coupling function is obtained as 

(6s - l)t 6s+3 



F(t) 



3£ ( 7 (6s - 1)kH 2 - 6s 2 (3s - 2)t 6s+v 



(3.6) 



According to Eq. (3.5) it should be A > 0, which is accomplished for s > 2/3 to 
get real scalar function <f>(t). In the special case of s = 2/3 (A = 0), corresponding 
to presureless matter dominance, the expression for the scalar field simplifies to <fi = 
2/37 1 / 2 t~ 3 / 2 , and the explicit dependence of the coupling function on the scalar field 
becomes 

_8_ /2T 
81^ 2 V 3 



F(<f>) 



2\ 1/3 



1 



03+1/3 



(3.7) 



8 



If we take the integration constant 7 = in Eq. (3.4), then in order to be consistent 



with t/> > 0, s should be in the interval 1/6 < s < 2/3. In this case integrating in Eq. 



(3.4) gives the following solution for the scalar field 

t 



+ 4>i log 



6s 2 (2 - 3s) 



(3.8) 



t ' ri V (6s-l)n 2 

with s in the interval 1/6 < s < 2/3. Note that the particular value s = 1/2 
corresponding to radiation dominance is an allowed solution. The coupling function 



in this case is obtained by making 7 = 0, s — 1/2 and replacing t from Eq. (3.8) in 
Eq. (3.6), giving the result (assuming O = at t = 1) 



K 



(3.9) 



In any case, the solution (3.4) is more appropriate for describing the cosmological 



dynamics at early times, when the radiation or matter content are dominating. 
B. Power-law solutions with potential V((j>) 7^ 



Multiplying the Eq. (2.17) by 0, defining the function ip for <fi 2 , and replacing the 
product F(<f))<f) 2 from Eq. (2.15), the Eq. (2.17) reduces to first order equation with 



respect to the variables ip, H and V, and can be written as 
H'j- + (qH 2 -H^jip + 2#^ - 2 (3H 2 + H ) V + 6^ f W 1 2H 



(3.10) 



Note that in Eq. (3.10) the functions 6(x) and V(x) are not coupled, and we can 



use this fact in order to simplify the Eq. (3.10), by imposing a restriction on one 



of these functions (in other words, they may satisfy separate equations by imposing 
appropriate restrictions on one of the functions). Using this fact, we can limit the 
model to the class of potentials that satisfy the restriction 



dV ( 9 • \ H 2 1 

H— — [3H + H)V + 3— 3H 2 + 2H 



, , . . , . _.. , (3.11) 

dt V J k \ ) v 

then, the equation for the field ip significantly simplifies, but still giving an interesting 

dynamics as can be seen bellow 



+[QH 2 - II ) r 







(3.12) 
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Let us assume the solution a(t) oc t p and replace in Eqs. (3.11. 3.12). Is is easy to 



check that a particular solution to Eq. (3.11) is given by 



V(t) 



3p 2 (3p - 2) 1 
k 2 (3p + 1) £2 



3p-l 



(3.13) 



where C is an integration constant, and the solution to Eq. (3.12) is 

-(6p+l) 



(3.14) 



where is the integration constant. The scalar field from Eq. (3.14) is given by 

-3p+l/2 



2y^ot 
6p — 1 



(3.15) 



with these solutions the scalar potential (3.13), and the coupling function from (2.15) 
can be written explicitly in terms of the scalar field as follows 

3p 2 (3p - 2) 



V{<j>) 



(3p + 1) kH 2 



4 

6p-l 



and 



F(<P)=[(^o) (3P+1)]- 1 (J 
which give rise to accelerated expansion for p > 1. 



2(6p+l) 

6p-a 



18^p 2 



4 

6p— 1 



(3.16) 



(3.17) 



4 Cosmological solutions with dynamical equation 
of state 



Here we study the full model with the restriction (2.14), and will use the freedom in 
choosing the coupling function in order to obtain a cosmological dynamics closer to 
the one expected from observations. 
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4.1 Scalar field dominance 



Below we study various known proposals for dark energy, including the cosmological 
constant and matter (named cold dark matter) ACDM, varying cosmological constant 
VCC and other models. First we consider the case of scalar field dominance (i.e. 
pm — 0). Defining the variables x for log a and 9 for <fi' 2 (""' denotes the derivative 



with respect to x), the Eq. (2.17) can be written as (after multiplying by 



~ (H 2 6) + 3H 2 9 + ^ + %H 2< ^F9 + 3£H 4 ^-(F9) + 18£H A F9 = (4.1) 

\JjJb \JLJU KajJu \JjJb 



From Eq. (2.15), changing to the variable x, we can write the product F(f)' 2 = F9 as 
following 

F9 = ^-JL-JL (4.2) 



3£k 2 H 2 18£H 2 9£H 4 



taking the derivative of Eq. (4.2) and replacing F9 and d{F9)/dx into Eq. (4.1), we 
arrive at the following equation involving 9, H and V 



A d9 „ / 9 dH 2 \ ~dV 
2H 4 — + H 2 12H 2 + — — 9 + AH 2 — 
dx \ dx J dx 



2[6H 2 + 



dH 2 

dx 



V 



H 2 

+ 12— ( 3H 2 



dH* 
dx 



(4.3) 







In this manner, we obtain a first order differential equation for the functions 9, H 
and V. In what follows we will introduce the scaled Hubble parameter H 2 = H 2 / Hq 
and the dimensionless scalar potential defined as V = k 2 V/Hq. Note that as in the 



above studied case, in Eq. (4.3) the functions 6{x) and V(x) are separated (in the 



sense that are not coupled), so we may proceed as in previous case by imposing a 
restriction on the scalar field potential V(x), given by the equation 



~~dV ( ~ 9 dH 2 , 
2H 2 — - \6H 2 + —— ) V + 6H 2 
dx \ dx 



3H 2 + 



dH 2 

dx 







(4.4) 



which simplifies the Eq. (4.3): 



2H' 



dx 



+ 12H 1 + 



dlP 

dx 



9 = 



(4.5) 
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where we have used the new defined scaled Hubble parameter H and potential V. 
ACDM solution 



In order to consistently solve the Eqs. (4.4) and (4.5), we propose the following 



expression for the scaled Hubble parameter H 2 



H z = Vt m oe 3x + Qa, ^mo > 0, > 



(4.6) 



where f2 m0 can be assumed to be the matter density parameter (this assumption is 
supported by constraints imposed by different observational datasets), and fl\ can 
be related with the cosmo logical constant. This is the known ACDM model (note 
that in this case the solution with cosmological constant is obtained as a result of 
the dynamics of the present model). It should be stressed the known form of the 



expansion factor, which from Eq. (|4.6|) is given by 



a(t) 



a- 



1/3 r 



sinh 



1 2/3 



(4.7) 



This solution describes a matter dominated universe at early times (a(t) oc t 2 ^ 3 ), and a 
cosmological constant dominated phase at late times (a(t) oc exp (3/2y/Q\Hoty). This 
model as is well known, is in excellent agreement with different sets of observations. 



The constants Q m o and Q A i n (4.6) satisfy the flatness condition at x = 



(4.i 



Replacing H 2 in (4.4), after integration we obtain the dimensionless scalar field po- 
tential 

-3x\ 1/2 



V{x) = 6^ A + 



6Q\e 3x 



n 



mO 



Ce 3x (Q A + Q m0 e- 3x y 



And replacing H 2 in (4.5) we get the following expression for 6 



6(x) = <p 



12 



-9x/2 



.1/2' 



(4.9) 



(4.10) 



(tt A e 3x + VL m0j 

where is the integration constant. Integrating the square root of this last equation, 
we obtain the scalar field as (considering the (— ) sign root) 



[x 



9a 



1/2 



rnO 



[Qa + ^moe 3x ) 



3^\3/4 



(4.11) 
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Finally, the coupling function F is found by replacing the Eqs. (4.6), (4.9) and (4.10) 



in (4.2). According to the expressions (4.11 ) and (4.9), in this case is possible to solve 



explicitly the potential in terms of the scalar field, giving the result 



V 



70 4 / 3 - n A 

where we assume the restriction 70 4//3 > 
VCC-type solution 



7 = Q 8/3 Oo 



4/3/3-2/3 



(4.12) 



We will consider the following solution, which generalizes the H power-law model 
in the context of the variable cosmological constant theories (VCC) [60] , |61j . 



H 2 = (Ae- ax + B)' 



(4.13) 



where A, a, B and /3 are constants that can be fixed by the expected behavior 
of the main cosmological parameters (The equation of state parameter (EoS) w, the 
deceleration parameter q, etc.), or can be constrained using the different observational 



datasets. Replacing (4.13) in (4.4), we obtain the scalar field potential of the form 



V(x) = Ce 3x (Ae- ax + Bf 2 + {) ' 



(a/3 + 6)(a/3-2a + 6) 



x 



W(a/3 + Q)e ax 2 F X 
(6a + 3a/3 - 2c? (3 + a 2 /3 2 - 18) 2 F X 



1 _3_/3 2 _3_/3 _Be« 
2' a 2' a 2 ' A 



~2'~a~2' ~a~2' A~ 



(4.14) 



where 2-F1 is the hypergeometric function. Integrating the Eq. (4.5) we obtain the 
expression for 9 

0{x) = 6 e {a P- 12)x/2 (A + Be ax )~ p/2 (4.15) 



where C and #0 are integration constants. After integrating the square root of (4.15), 
gives the scalar field 

49l /2 A-V* 



(j)(x) 



a/3 -12 



_ e (a0-12)*/4 2i?i 



_3 /3/3 3 ^ Be a 
a 4'4' a 4' A 



(4.16) 
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Finally, replacing H 2 , V and 9 from (4.13)-(4.15) in the generalized Friedmann equa- 



tion (4.2), we obtain the expression for the coupling function F, thus completing the 



solution to the Eqs. (4.2) and (4.3) under the restriction (4.4). Here the parameters 



A and B must satisfy the condition A + B = 1 at x = 0, as follows from (4.13). If 
we take A = f2 m0 , a = 3(1 — |) and (5 = 2/(2 — n) in (4.13), then we arrive at the 



solution to a variable cosmological model of the H power-law type with matter-dark 
energy interaction [60J. 



FRW with 2D topological defects 



An analytical expression for the potential and the scalar field can be found for the 
Hubble function of the form 

H 2 = n m0 e- 3x + n T0 e- x , (4.17) 

which describes the FRW cosmology with 2D topological defects, with the EoS of the 

formpx = — 2/3px 
for the potential 



. Replacing in Eq. (4.5 ) we obtain the following expression 

n T0 e 2 * 



V(x) = Ce^ 2 {Q m0 + Q T0 e 2x ) + ^to U 2 * + ^) 2ft 

V ^ 'mO / 



5 1 

'4' 2' 



1 

4' 



(4.18) 

and the function 9 is given by 

9{x) = V29 e~ 9x / 2 (n m0 + n T0 e 2x ) 1/2 (4.19) 

However, when a fixed value for Q m Q is considered, then while Q m o < 0.22 the model 
with 2D topological defects gives an appropriate description of the EoS, according to 
the current observations (i.e. this model favors a low density universe 



FRW with phantom dark energy 



The following expression describes the FRW solution with phantom dark energy, 
with EoS for the phantom sector of the form px = ^xpx [HI] 

H 2 = Q m0 (l + zf + Q P0 (1 + z )^x) 



(4.20) 
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where wx < — 1 and the parameters Q m o and Q po satisfy the restriction at the present 
epoch Q mQ + Qp = 1. Turning to the x- variable and replacing this proposal in Eq. 



(4.4), we get the potential 

V(x) =Ce 3x/ Wn m o + tt P0 e- 3 ^- 



3 + wx 



3(l+wx)i 




i F 1 



PO 



1 3 + wx —3 + wx 
2' 2wv 2w 



^mO e 



3wjf x 



X 



(4.21) 



where C is the integration constant. From Eq. (4.5) it follows the expression for the 
function 6(x) 



B{x) = e 



-9/2x 



(Q m0 + Q P0 e-^) 1/2 



(4.22) 



replacing H 2 , V and 6 from (4.20)-(4.22) in (4.2), we obtain the expression for the 
coupling function F, which allows to completely solve the generalized Friedmann 



Eq. (4.2) and equation of motion (4.3), under the restriction (4.4). In this case 
it is interesting to analyze the EoS parameter, which as expected presents quintom 
behavior [S3] . To this end, we write the effective equation of state parameter in terms 
of the redshift 



w 



eff 



1 (1 + z) d(H 2 
3 #2 dz 



(4.23) 



replacing H 2 from Eq. (4.20) we obtain 

w e// (z) = 



1 



1 + z) [3O m0 (l + z) 2 + 3(1 - n m0 )(l + wx)(l + z)^^- 1 ] 



3 [n m0 (i + zf + (i - n m0 )(i + z)3(i+wx)] 



(4.24) 

where we used flpo = 1 — fi m0 - Evaluating this expression at the present epoch [z = 0) 
gives w e //o = wx — fi m oWi, showing that this model crosses the phantom barrier and 
presents quintom behavior for Q m0 < (|wx| — 1)/ |wx|- In Fig. 1 we plot w e ff(z) for 
two representative values of Wx, taking Q m0 = 0.27. 
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1+z 

Fig. 1 The effective equation of state w e ff versus redshift, for two values ofwx ■ 
red-(wx = —1-2), blue-(wx = —1.5). The continue line w e // — — 1 corresponds to 
the cosmological constant. The dashed line (w e // = —1/3) intersects the curves at 
zt ~ 0.65 — 0.72. Note that the smaller the \wx\is, the earlier the redshift transition 

occurs 

4.2 Solutions with both scalar field and matter contributions 



Let us now focus in the solutions with the contribution of the matter term. Here we 



will consider the ACDM and a more general model of dark energy. From Eq. (2.15) 
and changing to the variable x, we can write the product H 2 F9 as follows 



H 2 F9 



1 



3^ 2 



9 



V 



Pmoe 



-3.T 



%H 2 %H 2 



(4.25) 



where we used the expression for the presureless matter component p r 



Pmoe 



-3x 



Taking the derivative with respect to x in Eq. (4.25), and replacing F9 and d{F9)/dx 



into Eq. (4.1), we arrive at the following equation involving 9, H and V 

dH 2 ~ 



a d9 o 
2H 4 — + H 2 
dx 



(l2H z + 

12 t9 ( 9 dH< 
+ —H 2 [3H 2 + — 
k \ dx 



dH 2 

dx 



9 + 4H' 



,dV 



PmO 3if + 



dx 
dH 2 



dx 



V 



dx 



(4.26) 



-3x 
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This is a first order differential equation for the functions 9, H and V. Starting from 



Eq. (4.26) we will consider different possibilities of cosmological scenarios. In order 



to integrate the equation (4.26), we use again the fact that the functions 9 and V 



are separated, and thus we may constrain the scalar field potential. It is important 

to note that we have different possibilities for the restrictions we can impose on the 

potential. 

First restriction 

Here we consider the following restriction on the potential, using the above defined 
H and V 



dx 



6H 2 + d ^)v + QH*(3H*+ dH2 



dx 



dx 



3Q 



mO 



3H 2 + 



dlP 

dx 



-3.c 



which after the replacement in (4.26), simplifies the equation for 9 : 



2H^ 

dx 



UH 2 + | 9 

dx 





(4.27) 
(4.28) 



where f2 m0 is the matter density parameter defined as f2 m0 = « 2 p m o/ '(3Hq). Note 



that this equation for 9 is the same as obtained in (4.5). 



ACDM solution 



Let us propose again the ACDM solution for H given by (4.6), and evaluate V 
and 9 from the Eqs. (4.27) and (4.28). Replacing H 2 in (4.27), after integration we 



obtain the following expression for the dimensionless scalar field potential 



V{x) = 3Q A + Ce 3x (fi A + Q m0 e- 3x ) 



3x\l/2 



(4.29) 



with C as an integration constant. From Eq. (|4.28|), it follows the expression for 9 

(4.30) 



9{x) = </>'• 



-9x/2 



(^ A e 3x + Q 



mO) 



,1/2' 



where is the integration constant. Integrating the square root of Eq. (4.29), we 



find the same solution for the scalar field obtained in (4.11 ). The explicit dependence 
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of the potential on the scalar field is obtained from Eqs. (4.11) and (4.29) 

cy/ 2 fi mO 2/3 



V = 3ft A + 



7< /> 4 /3_ft A • 



(4.31) 



From Eqs. (4.9) and (4.29) it follows that the dark energy potentials (without and 
with matter content) describing the ACDM model, increase as the Universe expands 
(x — > oo). 



Phenomenologically the desired solution 



In this case we consider an interesting solution for the Hubble parameter, which is 
usually introduced phenomenologically to describe the dark energy, and successfully 
fit the observations (see [60] ) 

H 2 = Q m0 e~ 3x + Q a e~ ax = ft m0 (1 + zf + ft Q (1 + zf 



where ft m o and ft Q are subject to the restriction 

ft m0 + ft Q = 1 



(4.32) 



(4.33) 



replacing the solution (4.32) in Eq. (4.4) we obtain the following expression for the 
scalar field potential 



V(x) = Ce-^ a - 3)x (n m0 e ax + n a e 3x ) 



3x W2 



3ft a (a — 3) 



ft m0 (2a + 3)(4a-3) V H 



1 + ^5L e (3-a)x , x 



'mO 



1/2 



2(2a + 3)ft a 2 F 1 



3-4a 3 9-6« ft 



(4a - 3)Q m0 e^ x 2 F, 



6 — 2a 2 6 — 2a ft m0 

2a + 3 3 3 -4a ft 



+ 



2(a-3)'2'6-2a' ft, 



mO 



and replacing (4.32) in Eq. (4.27) gives the solution for 9 

e(x) = 6 



which in turns gives the scalar field 

9 



x 



) = 0oe" 9 * /4 2 F 1 



Vtt a e 3 * + ft m0 e c 



1 + 



ft 



4(a-3)'4' 4(a-3)' ft 



«_ e (3-a)a: 
mO 



(4.34) 



(4.35) 



(4.36) 
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where 



40, 



1/2 



I) /(9^mO 



The parameters f2 m o and a in (4.32), can be constrained by demanding that the 



current value of the EoS parameter be Wq = — 1. To his end, we replace the scaled 



Hubble function (4.32) (given in terms of the redshift) in the expression for w e ff 
(4.23), and impose the condition w e //(0) = w 
relation 



1. After that, we obtain the 



3Q 



a 



mO 



^mO ~~ 1 

If we take for instance the value Q m Q = 0.27, then a > 
that wo = —1.1, in which case for Q m Q = 0.3, a takes the value a 
model presents quintom behavior. The behavior of the effective EoS parameter, for 
both cases is shown in Fig. 2 



(4.37) 

-1.11. We can also demand 
-1.71 and the 



^ -0.8 




0.6 0.8 1.0 



1.8 2.0 



1+Z 



Fig. 2 The effective equation of state w e ff versus redshift, for two sets of 
parameters: red-(fl m0 = 0.27, a = —1.11), blue-(fl m0 = 0.3, a = —1.71). The 
continue line w e // = — 1 corresponds to the cosmological constant. The dashed line 
( w e// = ~ 1/3) intersects the curves at zt ~ 0.57 — 0.68. Note that the increase in 
Q m o delay the onset of the redshift transition. 

Note that the curves in Fig. 1 and 2 correctly describe the matter dominance at early 
times, the transition to accelerated expansion at z ~ 0.57 — 0.72 (see [HI]), and the 
current accelerated expansion with the equation of state very close to —1. 
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Concerning the time variation of the gravitational coupling, and applied to the solu- 
tions found here, we can proceed in the same way as in [23] and [21]. The effective 



gravitational coupling from (2.15) can be written as 

G 



G 



eff 



(4.38) 



1 - 3£k 2 FH 2 9 

H 2 d>' 2 = H 2 9. Then, in terms of x the time 



where we used k 2 = 8tcG and <fi 2 
variation of the gravitational coupling can be written as (taking into account that 
d/dt = Hd/dx) 

Geff 3^ 2 



d 



1 - 3£k 2 FH 2 6 dx 



(FH 2 6)H. 



(4.39) 



G eff 

Replacing the product FH 2 9 from Eqs. (4.2) or (4.25), and evaluating at the present 



time (as = 0), the Eq. (4.39) can be written as 

Geff 3/(C, 9 ) 



G 



eff 



i=0 



1-3<7(C,0 C 



-H 



(4.40) 



where the main parameters of the model appearing in H 2 have been fixed, and the 



only free parameters present in (4.40) are the integrations constants C and 9 . This 
defines the functions / and g through the expressions f(C,8) = £k 2 FH 2 9 and g = 
£k 2 -£;(FH 2 9), valuated at x = and at the values of the parameters that have been 
previously fixed. We can meet the constraints on the current value and the time 
variation of the gravitational coupling |65j , by appropriately defining or constraining 
the constants C and 9 . Here we illustrate the important case of the ACDM solution 



( |46) with ( [4729] ) and gTjjjp . Replacing ( [46) , g2g and ( [4730] ) in Eq. g25j, and 
then in (4.39) (making the needed changes to H 2 and V in (4.39)), and evaluating at 
the present time, it is obtained 



G 



eff 



G 



eff 



5 Jl 



K 2 9 n 



2 = 



1 + 



1 + 



10 



Hn 



(4.41) 



3(2+n m0 )_ 

In this case we can satisfy the requirements of the time variation of the gravitational 
coupling, by defining the constant C as C 



5k 2 9 



annulling the expression (4.41). 
And by choosing 9$ = Hq, we make the second term in the denominator (and there- 



3(2+O m0 ) 



fore in the denominator of (4.38) which affects the current value of the Newtonian 
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coupling) extremely small, all in complete agreement with the observed constraints. 



Other restrictions on the potential 



Besides the above restriction (4.27), and due to the extra term in Eq. (4.26) 



related to matter, we can consider other possible restrictions on the potential V(4>). 



Consider for example the same restriction given by Eq. (4.4). After the replacement 



in (4.26), the corresponding equation for the function 6(x) becomes 

dH 2 



2H^ + H 2 

dx 



12H Z + 



dx 



9-3tt 



3H Z + 



dJP 

dx 



-3x 



(4.42) 



Thus, we keep the same solutions for the Hubble parameter and the potential as 



obtained in subsection 4.1, but according to (4.42) the equation for the scalar field 



becomes more complicated and some solutions can not be expressed analytically, and 
should be solved numerically. 

Other possible restrictions on the potential are given by 



or 



dV { ~ 2 dH 2 

2H 2 - \QH 2 + — 

dx \ dx 



2H^ 

dx 



V - 3n m0 3H 2 + 



dH 2 

dx 



-3x 







QH 2 + 



dH 2 

dx 



V = 0, 



(4.43) 



(4.44) 



Hence, starting from the model with the action (2.1 ) we obtained a variety of solutions 



which not only give the accelerated expansion, but also give the correct transition from 
decelerating phase to accelerating phase (for transition deceleration-acceleration in 
f(G)-gravity see [66]), in agreement with the observational data. The fact that the 



present model contains the type of solutions as described above in subsections 4.1 



and |4.2[ indicates that the scalar field with non-minimally kinetic coupled terms, as 
presented here, could explain the dynamical origin of these (mostly) phenomenologi- 
cal models of dark energy. On the other hand, in order to obtain the above solutions 



see section 4.1) it was not necessary to introduce the matter term in the action. The 



introduction of matter term translates into extra term in the Eq. (4.26), allowing for 
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different restrictions on the potential, which may enrich the dynamical possibilities 
of the present model. 



5 Discussion and conclusions 

We have considered a model of scalar field with kinetic terms coupled non-minimally 
to the scalar field and to the curvature. We found power-law solutions without and 
with a potential, giving rise to different power-law scenarios which include radiation 
dominance, matter dominance and also describing accelerated expansion. Due to the 
complexity of the system of equations which describes the dynamics of the model, and 
thanks to the presence of the coupling function F ((/>), we could impose a restriction 



on the potential through Eq. (4.4) or (4.27) (if explicitly matter content is added), 



which allowed us to find a class of solutions that lead to dynamical effective EoS. We 
proposed different solutions for the Hubble function, corresponding to known impor- 
tant models of dark energy, like ACDM, variable cosmological constant model VCC, 
FRW with two dimensional topological defects, FRW with phantom dark energy, and 
so on. These solutions properly describe the accelerated expansion and the transi- 
tion from the decelerated to accelerated phase, according to the current astrophysical 
observations, and also account for the cosmic coincidence (i.e. that guarantee the 
matter dominance at early times (high redshift), and the dark energy dominance at 
late times). This means that the present model could explain the dynamical origin 
of many phenomenologically viable dark energy models. 

One relevant fact of this model, is that it contains solutions with EoS that cross the 
phantom frontier, showing that the present model with non-minimal kinetic coupling 
is able to describe the quintom behavior. We also illustrated the time variation of the 
gravitational coupling for the relevant case of the ACDM solution. But in general we 



can proceed with the Eq. (4.39), and in order to satisfy the current restrictions on 
the Newtonian coupling, we can use the freedom in the integrations constants C and 
#o, in order to control the the actual value of the Newton's constant and its time vari- 
ation. These conditions can be satisfied by imposing the inequalities f(C, 9 ) < 1CT 1 
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and g(C,9 ) < lO" 5 . 

In conclusion, the derivative couplings provide an effective scalar density and pressure 
which might play an important role in the explanation of the dark energy or cosmo- 
logical constant. The present model shows a wide variety of acceptable solutions, and 
particularly we found solutions for the Hubble parameter, used to describe the actual 
phenomenology of the dark energy, according to the current observational data. 
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